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Abstract. - 

We study the statistics of turbulent velocity fluctuations in the neighbourhood of a strong large 
scale vortex at very large Reynolds number. At each distance from the vortex core, we observe 
that the velocity spectrum has a power law "inertial range" of scales and that intermittency - 
defined as the variation of the probability density function (PDF) of velocity increments as the 
length of the increment is varied - is also present. We show that the spectrum scaling exponents 
and intermittency characteristics vary with the distance to the vortex. They are also influenced 
by the large scale dynamics of the vortex. 
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' Introduction. - Much efforts have been devoted to the study of high Reynolds number 



turbulence, assuming the properties of local homogeneity and isotropy. Under these assump- 
tions it has been shown jl| that, in between the integral scale at which energy is fed into the 
flow and the dissipative scale at which viscosity smoothes out the velocity gradients, there 
exists an "inertial range" where: (i) the velocity spatial power spectrum follows a power law 
E{k) ~ e 2 / 3 £;~ 5 / 3 , (ii) the energy transfer rate e is related to the third order structure function 
Ss{r) = (5u 3 (r)) — — fer. The first point has been verified experimentally in numerous 
experiments and constitutes the main success of Kolmogorov's K41 mean field theory ||. 
The second point is generally admitted although experimental verifications are rarer In 
addition, it has been observed that the probability distribution of the amplitudes of velocity 
increments depends on the width of the increment. It varies from a gaussian PDF at integral 
scale to the development of stretched exponential tails at the smallest scales. This feature 
defines intermittency, and traces back to the non-uniformity in space of the energy transfer 
rate. These three observations are the building blocks of the studies of homogeneous, isotropic 
tubulence (HIT). 

It is the aim of our work to study experimentally how these observations are modified in a 
situation of inhomogeneous turbulence. Partial studies have been made near weak vortices [Q, 
in inhomogeneous wakes (|, || or near a boundary layer j^, ||. Here we study turbulence 
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Fig. 1. - Experimental setup. The fluid is confined in a cylindrical vessel of height 40 cm and 
radius 11.7 cm. The disks, of radius R = 9.8 cm are set a distance H = 30 cm apart. They are 
driven by two d.c. motors whose rotation rates Oi , O2 are adjustable in the range [15,45] Hz. 
The hot-film probe is located just above the midplane, at variable distance r from the axis of 
rotation. 



in the vicinity of a very intense large scale vortex. Our motivations are twofold: first, most 
high Reynolds numbers flows are inhomogenous and the question of the "universality" of the 
features of HIT should be addressed. A second motivation is that it has been proposed that 
such coherent structures play a significant role in the statistical properties of turbulence [|[ 
H- 

Experimental set-up and measurements. - The experimental apparatus and measurement 
techniques are those described in fl3|| . A strong axial vortex is formed in the gap between two 
corotating disks - fig. 1. The control parameters of the flow are the disks angular velocities: 
fii, and the mean integral Renolds number: Re = R 2 ^/Q\ + £l\/v. Essentially two regimes 
are observed. The first one, hereafter labeled (GR), is when Q,\ ~ f^: the core of the flow 
has a strong global rotation and a stable large scale axial vortex is present. On the other 
hand, when one disk rotates at a much faster rate than the other, the flow is dominated by 
differential rotation; a strong axial flow is produced and one observes intermittent sequences 
of formation and bursting of large scale axial vortices - in the following, we label (DR) 
this regime. In each case, the integral Reynolds number of the flow exceeds 5 .10 5 so that the 
flow is 'turbulent'. 

Local velocity measurements are performed as in [fL3|. The hot-film probe is set parallel to 
the axis of rotation, so that it measures u r g — \Ju 2 + Ug (here, (r, 8, z) are the conventional 
cylindrical coordinates with z parallel to the rotation axis). The sampling frequency Fs is 
78125 Hz for measurements in the (DR) regime and 39062.5 Hz in (GR) - In each case, at 
least 10 7 data points are collected. We stress that although the flow has in most locations a 
well defined mean value u of the velocity u r g and fluctuation levels it rms comparable or less to 
thoscs reported in jet turbulence, we do not attempt to use and justify a 'Taylor hypothesis' 
here and we analyse our result as time series. 

Spectra. - A first question regarding such a turbulence concerns the range of scales of 
motion and the energy distribution among them. In fig. 2 we show the power spectra of 
velocity time series recorded at different distances from the rotation axis. It can be observed, 
as a first sign of the flow inhomogeneity, that the spectra depend very much on the location 
of the probe and on the flow configuration. 

In all cases, we observe the existence of an 'inertial range' u E [u>i, ujt] where the spectrum 
follows a power law E(oj) oc u>~ a . However both the range u>i/u>t and the exponent a vary. 
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- In (GR), a increases from about 1.65 in the outer region, to 2.50 at the rotation axis; 
meanwhile the scaling domain is also reduced. A more detailed analysis of the spectra shows 
that the energy in the larges scales remains approximately constant while both the frequency 
and the energy at the end of the 'inertial range' decrease with r. As a result, the proximity 
of the stable vortex structure reduces the range of scales of the turbulent motion and their 
energy content. This is consistent with earlier findings |TJ, [l^] which have shown that global 
rotation tends to decrease the intensity of small scales motion. 

- In (DR), a decreases from again about 1.65 in the outer region, to 1.1 at the rotation axis. 
The scaling range does not vary. Here, one observes that the energy in the dissipative region 
is almost independent of the position r of measurement. It is the energy measured in the large 
scales that decreases with r. 

We thus find that the dynamics of the large scales influences the entire range of scales. It is 
surprising to observe that some self-similarity is retained since the spectra follow a power law 
(most curves exhibit almost a decade of scaling behaviour) : there is no scale separation where 
one would observe the pathologies of the injection at large scales and where the turbulence 
would recover the characteristics of HIT at small scales. 

Third order structure function. - The energy transfer across the scales is usually illustrated 
by the behaviour of the third order structure function. In HIT, this is justified by the Karman- 
Howarth relationship. But an anistropic equivalent can be established Q. We show in fig. 3 the 
evolution of the third order moment of the velocity increments 63(7") = {(u(t) — u(t + t)) 3 )^ 
computed at different distances from the rotation axis. 

- In (GR), we observe that in the outer region, S3 (r) is a negative, bell-shaped curve. Negative 
values of the skewness indicate that the cascade proceeds in the forward direction, from the 
large scales to the dissipative ones, as in HIT. At such limited Reynolds numbers, one cannot 
expect a plateau in plots of Ss(t)/t (even less with value —4/5), but the maximum of \S3(t)/t\ 
is consistent with the value of the power consumption per unit mass e, measured from the 
electric consumption of the motors p6| . On the contrary going towards the axis the amplitude 
of S3 (t) diminishes rapidly. Near the vortex core, the skewness eventually becomes positive 
in an intermediate range of increments. This indicates that the cascade to small scales is first 
reduced, and then reversed by the influence of the coherent rotation of the vortex. Energy is 
accumulated in the large scale where it stabilizes the dynamics of the large scale vortex. 

- This picture is quite different in the irregular (DR) regime. There, except at the rotation 
axis, the skewness has the behaviour traditionnally observed in HIT. The cascade proceeds 
forward; the maximum shows little variation with the location of the probe and has a value 
which is, again, consistent with the motors power consumption. At the axis, Ss(r) drops 
suddenly; even its sign becomes unclear and it may even be slightly positive in a range of 
scales. 

Altogether, the analysis of the third order moment shows that the structure and dynamics 
of the flow at large scale strongly influences the entire energy transfer across the scales. 

Intermittency. - A first observation is that intermittency is always present, regardless 
of the flow regime and probe location. Indeed, for every measurement, the PDFs of velocity 
increment change from a quasi-gaussian shape at integral scale to the development of stretched 
exponential tails at small scales, see fig. 4 (the PDFs shown in this figure are for signed 
quantities). 

However, as it has now become customary, we analyse here intermittency through the 
behaviour of moments of the absolute values of the velocity increments: 

S p (t) =< \5u(t)\ p >=< \u(t + r) — u(t)\ p > t . A feature common to all flow configurations 
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and every distance to the axis, is that the evolution of S p (t) is well modeled by a self-similar 
multiplicative cascade in the sense that one can write: 

log(5 p (r)) = H(p)n(r) , 

in a large interval of increments ][l7| ^8|. Here, H(p) is the Laplace transform of the cascade 
propagator Jisj] which reflects the hierarchy of amplitudes of the velocity increments and n(r) 
describes the speed at which the cascade proceeds along the scales. Experimentally, we have 
checked the above relation in the following manner: 

(i) the plot of one moment against an other yields an extended scaling region, with exponent 
H(p)/H(q) if one plots S p vs. S q - this is the ESS property ^0|. Since the decomposition of 
a moment S p (r) into the two functions H and n is defined up to an arbitrary multiplicative 
constant, we set H (3) = 1. Unlike the context of HIT, where this choice relies on the expected 
scaling of the third order structure function, it is here a priori arbitrary. 

(ii) once that first step is done, and the values H{p) obtained, one computes for every order 
the function n p (r) = log(S p (r))/H(p). We so verify that n p (r) is independent of p and obtain 
a unique function n{j) for every measurement. 

Following this procedure, the results are: 

- Up to experimental errors, the functions H (p) do not depend on the measurement location 
nor on the flow regime. The values (H(p) = 0.42, 0.70, 1.27, 1.51, 1.73 for p = 1, 2, 4, 5, 6) are 
those reported in HIT in numerous experiments || pl| . 

- The functions n(r) do depend on the measurement location and flow regime. They are shown 
in figure 5 (the curves are normalized so that n = at integral time delays). We do not discuss 
hereafter the functional form of the n(r) curves but the relative behaviour of the curves in 
each regime, when r is varied. In (GR) one observes that in the outer region a larger range 
of scales is covered with a given number of cascade steps n than in the vicinity of the vortex 
core. The slopes dn(r)/dr increase as r decreases, meaning that increasing number of steps are 
needed to cascade between any two given scales. Again, this is consistent with the idea that 
rotation prevents the 3D cascade and with our previous observation that it reduces the energy 
transfer to smaller scales. An inverse effect is observed in the (DR) regime: as measurements 
are performed closer to the rotation axis, ti(t) increases and is less steep, meaning that a 
wider range of scales is reached in a given number of cascade steps. This is consistent with a 
reduced slope of the velocity power spectra since such an efficient cascade provides an enhanced 
distribution of energy in the smaller scales. 



Concluding remarks. - Our foremost original observation is that some major caracteristics 
of HIT are retained in the vicinity of large scale coherent vortex structures: (i) the spectra 
develop an "inertial range" in the sense that power law scaling regions exist, (ii) the 3rd order 
structure functions yield consistent descriptions of the energy transfers and (Hi) intcrmittency 
is always present and can be adequately described by a self-similar multiplicative cascade 
model. Second is that the structure and dynamics of the large scales influence the entire range 
of scales of motion. The combined study of power spectra, third order structure function 
and cascades parameters H and n, give valuable informations that show that the turbulent 
cascade is reduced in the presence of global rotation and enhanced in the vicinity of very 
unstable vortex structures. Note that in the (GR) regime, eventhough the energy cascade is 
reversed and the spectra evolve towards a "—3" slope in the inertial zone, intcrmittency is 
observed: this a sharp difference with 2D turbulence |22 where no intermittency is detected 



although deviations from Gaussianity may be present [23]. 
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Fig. 2. - Power spectra of velocity times series recorded at distances (0, 0.5, 1, 1.5, 2, 2.5, 
3, 3.5, 4, 4.5, 5.5) cm from the rotation axis, (a): Q\ = O2 = 30 Hz; (b) : Oi = 40 Hz, 
2 = 12 Hz. 




Fig. 3. - Evolution with increment length of the third order moment of the veloc- 
ity, (a): (GR) regime; (b): (DR) regime. Ss(t) is computed at distances r = 
0.5(o), 1.5(0), 2.5(D), 3.5(A), 4.5(*) cm from the axis. The horizontal axis is in units 
of the sampling frequency. 



c. simand et al: inhomogeneous turbulence. 



7 




Fig. 4. - Evolution with increment length of PDFs of the velocity increments, (a): (GR) 
regime; (b): (DR) regime. The PDFs are computed at distance r — 2.5 cm from the axis 
and for increments r = (1, 5, 21, 49, 116, 271, 635, 1485)/Fs. Values of the velociy have been 
normalized by the rms fluctuation level. The curves are translated for clarity. 




Fig. 5. - Evolution of the cascade depth, (a): (GR) regime; (b): (DR) regime. n(r) is computed 
at distances r = 0.5(o), 1.5(0), 2.5(D), 3.5(A), 4.5(*) cm from the axis. 



